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. spheric neutrino data, especially at energies higher than approximately 50 GeV, 

where the effect of the tau mass is small. For example, a 100 GeV NUFACT can 
probe FC neutrino interactions at the level of few |e| < few x 10~^ at 99 % C.L. 
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standard neutrino-matter interactions (NSI) . We find that the sensitivity to flavour- 
changing (FC) NSI can be substantially improved with respect to present atmo- 
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1 Introduction 



A long baseline neutrino factory [1-3] offers a unique tool for addressing basic 
questions in weak interaction and flavour physics. One outstanding example 
is the quest for neutrino mass and oscillations, which touches fundamental 
issues related to Grand Unified Theories. Motivated by the great discover- 
ies of underground experiments [4-8] neutrino mass and oscillation searches 
have become the center of attention in particle physics research. Apart from 
being motivated on basic theoretical grounds [9-12] neutrino masses and os- 
cillations offer the simplest and most obvious way to account for the observed 
anomahes [13]. Nevertheless other mechanisms, based on flavour changing non- 
standard neutrino interactions have been suggested in connection with both 
solar [14-16] and atmospheric anomahes [17-20] as well as other astrophysics 
applications [21,22]. They can either provide alternative solutions [23] or else 
be severely tested by the data, in the atmospheric case [20]. They may arise 
in a number of theories beyond the Standard Model [24-27] , in particular, in 
most (but not all) models of neutrino masses [10]. 

Using neutrinos from an accelerator in order to obtain an independent con- 
firmation of the non-accelerator physics results, is of fundamental interest, as 
it will bring more light upon the issue of neutrino masses and oscillations. 
This has been the focus of a number of dedicated recent NUFACT studies [1- 
3,28]. Following the recent suggestion in [29] we propose the use of a generic 
neutrino factory (NUFACT) to probe non-standard neutrino-matter inter- 
actions (NSI). We show how indeed such an ideal NUFACT can improve our 
present knowledge of non-standard FC neutrino interactions well beyond what 
is presently attainable on the basis of the latest atmospheric results and dis- 
cuss the corresponding energy, luminosity, energy resolution and tau detection 
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requirements. We find that, for example, a 100 GeV NUFACT can probe FC 
neutrino interactions at the level of few < few~^ at 99 % C.L. without 
any assumption about tau charge identification. In contrast no improvement 
is expected on non-universal (NU) neutrino interactions beyond the present 
achieved sensitivity. 

In order to compare the NUFACT sensitivities to NSI with present atmo- 
spheric sensitivities, we will adopt the same approximation as in ref. [20], i.e. 
we neglect the possible NSI in the production and detection process of neutri- 
nos. It is well understood that NSI can be probed in a near detector with high 
accuracy (see e.g. [30]). However, the event rates in a near detector depend 
quadratically on the strength of the NSI, whereas exploiting the non-standard 
matter effects we obtain a linear dependence of the rates in a far detector. 
A combined treatment of NSI in production, propagation and detection would 
lead to intriguing interference effects and is beyond the scope of this letter. 



2 Interplay of Neutrino Oscillations and non-standard Interactions 



The Standard Model can be extended to add neutrino masses in a variety of 
ways [12]. In any massive neutrino gauge theory the charged current (CC) 
weak interaction is characterized by the lepton mixing matrix Kaj- This neu- 
trino mixing matrix arises from the unitary matrix (U) diagonalizing the neu- 
trino mass matrix and the corresponding unitary matrix (Q) diagonalizing the 
left-handed charged leptons {K — QC/) and can be written in the following 
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parameterization [10] 



C12C13 



S12C13 



-iScp 



-S12C23 - Ci2S23Sl3e*''cP C12C23 - Si2S23Sl3e'^^^ S23C13 ' ^ > 
■5l2'523 ~ Ci2C23Si3e*'^'^P " C12S23 — Si2C23Si3e*'^'^P C23C13 ^ 

where we see explicitly the usual three neutrino mixing angles Ou, O23, 6*13 and 
one CP phase 6cp- This is the analogous of the CP phase found in the quark 
sector, as the other two Majorana phases were set to zero, since they are not 
observed in standard total-lepton-number-conserving Q oscillations. 

The above 3x3 form applies if there are no SU(3)xSU(2)xU(l) singlet lep- 
tons, such as the simplest models where neutrino masses arise radiatively [32,33] 
In seesaw type schemes [9,11] the matrix K^j is rectangular and contains in 
general many more parameters: twelve mixing angles and twelve CP phases 
in the three generation seesaw scheme [10]. We assume, however, that singlets 
are all super-heavy so that the K^j matrix can be well approximated by a 
unitary 3x3 matrix and parameterized as eq. (1). This is in fact in agreement 
with the scale of neutrino mass indicated by present neutrino anomalies. 

All present neutrino data Q can be accounted for by eq. (1). The two mass 
splittings Ami, = Am^g and Am^j.^^^ = Am23 ~ Am^g as well as the three 
neutrino mixing angles are all determined by global fits of neutrino data [13] 
which indicate that two of the angles are large, ^13 being small due mainly 
to reactor results [36]. The recent SNO CC data [37] adds support for the 

^ They could be seen only in AL = 2 processes, such as discussed in [31]. 

^ Except for the LSND anomaly, which requires a light sterile neutrino. For recent 

discussions see [34,35] 
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so-called LMA solar neutrino solution [38], which previously came only from 
detailed solar recoil electron spectra [39]. Moreover, LMA is also consistent 
with the the observed SN 1987A neutrino signal [40]. Thus in what follows we 
will take the parameters appropriate to this solution. However the details of 
the solar neutrino oscillation parameters do not significantly affect our results. 

Many theories beyond the minimal SU(3)xSU(2)xU(l) model also lead to 
non-standard neutrino interactions. These include most models of generating 
neutrino masses, which are generically accompanied by NSI, such as the sim- 
plest seesaw type schemes [10,12], super-gravity SO(IO) unified theories [25], 
models of low energy super-symmetry with broken R parity [26] as well as some 
radiative models of neutrino mass [27]. Exceptional examples exist of situa- 
tions where FC interactions are unaccompanied by neutrino masses. Models 
involve neutral heavy leptons at weak scale [41,42] and some super-gravity 
SU(5) models [24]. Such non-standard interactions [10,14-16] can be either 
flavour-changing (FC) or non-universal (NU). 

In Refs. [17-19] the atmospheric neutrino data have been analyzed in terms 
of a pure — > Ut conversion in matter due to NSI. The disappearance of 
from the atmospheric neutrino flux is due to interactions with matter which 
change the flavour of neutrinos. A complete analysis of the 79 kton-yr Super- 
Kamiokande data, including both the low-energy contained events as well 
as the higher energy stopping and through-going muon events from Super- 
Kamiokande and MACRO was given in Ref. [20]. 

We therefore study an extended mechanism of neutrino propagation which 
combines both oscillation (OSC) and non-standard neutrino-matter interac- 
tions (NSI). In order to discuss the sensitivity of NUFACT to non-standard 
neutrino interactions we adopt the general Hamiltonian given by 
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containing both non-universal (NU) and flavour-changing (FC) interactions 
characterized by diagonal and off-diagonal entries in eq. (2). 

Note that, as usual, the matter potentials for neutrinos and anti-neutrinos 
differ in sign. In contrast, we assume that the new interactions are CP con- 
serving. As a result the epsilon values have the same sign for neutrinos and 
anti-neutrinos. In the most general case the non-standard interactions might 
violate CP and the resulting phases could therefore affect the evolution. We 
will not consider this more complicated case in the following discussion. 

In addition to the flve standard parameters (three angles and two mass split- 
tings, if CP conservation is assumed) which describe the oscillation among 
three neutrinos there are, in the present scheme, also the and e^^ para- 
meters characterizing the NSI of the neutrinos. Of the three possible channels, 
we choose to analyze in detail here only the — v^. transitions closely related 
to the atmospheric anomaly^. The others will be discussed elsewhere. 

The relative importance of masses and NSI in the propagation of neutrinos 
is difficult to predict from basic principles and it is rather model-dependent. 
From a phenomenological point of view, however, atmospheric data imply that 
NSI can only play a sub-leading role [20] in z/^ — transitions. 



^ For this reason we have neglected the e^^ and e^^ involving the first generation 
in Eq. (2) 
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In order to gain some insight in the interplay between oscillation and NSI it 
is useful to reduce the problem to a two neutrino case by taking the limit 
Am^2 ~^ 0- this case the rotation in the 12-subspace also drops out, and 
therefore also the CP- phase becomes irrelevant [10]. This approximation is 
quite accurate for the v^—Vr transition at the energies and baselines considered 
for a neutrino factory experiment. In this limit only five parameters remain: 
three OSC parameters $23, ^13 and Am|i and our two NSI parameters e^^^ and 
e' . Neglecting ^13 the transition probability is given by [29]: 



P(-. - -r) = sin^ , 



2E 



B = Ai3sm2e23 + 2e^rVf, 
C = Ai3 cos 2^23 + e;.^^/. 



(3) 



(4) 
(5) 



3 Simulating neutrino factory long baseline experiments 

In testing the effect of non-standard interactions in the i/^ — i/j- transition it 
is essential to have a detector which is able to identify u^- events with a high 
efficiency. We are fully aware that this is a very difficult goal to achieve in 
the design of a large detector (m ~ 10 kt). However, this work is partially 
intended to show the benefits of such a detector in probing physics beyond 
the standard model. We will assume a detector with a mass of 10 kt which is 
able to detect and identify u-,- interactions above a threshold of 4 GeV with a 
constant efficiency of 77 = 0.33. Basically there are the following observables: 
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Ve Vr n+{Ur) 
U/^ Ur n+{Ur) 



(6) 



As will be clear in section 4 the ability to identify the charge of the tau is not 
necessary for this particular transition and therefore not assumed. This leaves 
us with only two observables: 



(7) 



where n and denote the event numbers arising from the neutrino factory 
operating in the two polarities. 



In calculating the event rate spectra in a neutrino factory experiment and for 
the treatment of the matter profile we follow the description given in ref . [43] . 
For the ly-r appearance channel we use the cross-section given in [44], the 
cross-section is assumed to be one half of this. We will show that neglecting 
the tau mass [29] is not a good approximation especially for neutrino energies 
below 20 GeV. We also take the energy resolution of the detector into account 
by modeling it as a Gaussian, as described in ref. [45]. The neutrino factory 
delivers 2 • 10^° useful muon decays of each polarity per year for a period of 
5 years. The energy of the neutrino factory is indicated in each figure since it 
plays a crucial role in probing non-standard interactions. 



We now describe the Statistical Method we employ. In order to estimate the 
sensitivity to new physics we adopt the following definition of [45] 

= 2 (n+ - n+sc) + 2n+sc In ^ + 2 (n" - n5sc) + 2n5sc In ^ . (8) 



8 



where riosc stands for the event rates which are expected in the absence of 
NSI. This is readily obtained by leaving all parameters as in the calculation of 
except that e^^ and e'^^ are set to zero. Thus has two degrees of freedom, 
therefore a value of 9.2 corresponds to 99% CL . Considering only total rates, 
the sum over the energy bins is performed before is calculated, whilst for 
an energy spectrum this sum is performed after calculating for each bin. 

The above is suitable to investigate the possible sensitivity to the new 
effects arising from non-standard interactions. In order to get reliable sensi- 
tivity limits it would be necessary in general to take into account possible 
parameter correlations and to evaluate the Vt -appearance together with i/^- 
disappearance and f^-appearance. However for the — transition this 
complication is less relevant to the extent the relevant parameters sin^ 2^23 
and Am|]^ that could be correlated with the NSI parameters are already well 
determined by present atmospheric data. As a consequence our results are 
basically unaffected by taking into account these correlations. We have in fact 
verified this by an explicit statistical analysis similar to that in ref. [45]. The 
situation would be different for the v^, — transition since this mode is con- 
trolled by sin^2^i3 which is subject to much higher uncertainties. For this 
reason this mode will be discussed elsewhere [46]. 

4 Results 

In order to highlight the effect of the non-standard interactions, we show in 
figure 1 the change in the ratio of events (which arise from to v^. events 
(which arise from for different values of the FC parameter e^,-. If no NSI 
interactions are present this ratio is basically constant 0.5 (black solid lines) 
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since this transition is only very weakly influenced by ordinary matter effects 
due coherent forward scattering off the electrons [15]. The value of 0.5 simply 
reflects the ratio of the cross-section for v and v. The grey shaded bands show 
the Gaussian Icr, 2(7 and 3cr errors on the standard ratio in the absence of NSI 
neutrino interactions. The dashed lines indicate the deviation from this for 
different values of the FC parameter e^T. One sees that for at the per cent 
level the difference is rather significant as long as the basehne is shorter than 
about 7000 km. For the left hand panel with an muon energy of 20 GeV these 
errors increase drastically with the baseline. This is due to the geometrical 
loss of flux at large distances. Comparing the two panels one can easily see 
the importance of high energies in order to obtain optimal sensitivity to the 
new physics. We have fixed the OSC parameters as follows: sin^ 2^^i2 = 0.78 
and Am|i = 3.3 • 10~^ eV^ (suitable to account for the LMA solution of the 
solar neutrino anomaly), sin^2^23 — 0.97 and Amg^ = 3.1 • 10~^eV^ (suitable 
to account for the atmospheric anomaly) and sin^ 2^^i3 = 0.02 in agreement 
with reactor results. Finally we have assumed CP conservation, bcp — ^ and 
also exact universality, e^^ = 0. The measurement of this ratio would require 
charge identification of the tau. It is shown here only for illustrative purposes. 

Note that the ratio in Fig. 1 contains only one part of the information con- 
tained in the event rates. For this reason we will use the as defined in 
equation 8 in order to calculate the sensitivity bounds to non-standard inter- 
actions. 

Before we do that let us highlight the important role played by present at- 
mospheric data by presenting Fig. 2. For the coming plots we use a base- 
hne of 732 km and a muon energy of 50 GeV. All other parameters are kept 
fixed as previously to: sin^ 26*13 = 0.02, sin^ 26*12 = 0.78, ^cp = 0, Am|i = 



10 



= 20 GeV E/, = 50 GeV 




10^ 10' 10" 10^ 10' 10* 

L [km] L [km] 



Figure 1. These figures show the ratio of observed Vr events from to the observed 
Vt events from as a function of the basehne. The grey shaded bands indicate the 
Gaussian la, 2a and 3cr statistical error on this ratio. The black solid line indicates 
the OSC prediction whereas the dashed lines indicate the deviation from this for 
different values of the FC parameter e^,-. The other parameters are sin^ 2^12 = 0.78, 
Am^i = 3.3 • 10-5 eV^, sin^ 26*23 = 0.97, Am^i = 3.1 • lO'^ eV^ sin^ 26*13 = 0.02, 
5cp = and e^.^ = 0. 

3.1 • 10-^eV^ Amii = 3.3 • IQ-^eV^ and e^^ = 0. In the left hand panel the 
dependence of the event rates for the — transition (solid line) and for 
the — Vt transition (dashed line) on the FC parameter e^T is shown for a 
fixed value of sin^ 2^23 = 0.9. For very small e^,- values the event rates in both 
channels are nearly independent of e^r and their ratio simply reflects the ratio 
of the cross sections. For increasing values of e^,- we now see that neutrinos 
and anti-neutrinos behave in a opposite way. This due to the fact that Vf has 
a different sign for neutrinos and anti-neutrinos. This behavior is also what 
would be expected from a linearized version of eq. 3 as given in ref . [29] . At 
e^T ~ 0.01 however this simple picture breaks down, the non-linearities of eq. 3 
become very important. Note that the transition probability only depends on 
B^. If the two terms contributing to B become of the same order of magni- 
tude, i.,e. Ai3 sin 2^23 — 2e^i-V/ then the difference between the sum of the 
two and their difference becomes maximal. Therefore the ratio of the anti- 
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neutrino rates to the neutrino rates becomes minimal. Increasing e^r further 
makes the difference between neutrinos and anti-neutrinos smaller again. For 
large values of e^^ the oscillation term A13 sin 26*23 becomes negligible and e^j- 
plays both the role of a mixing angle and provides the leading contribution 
to the mass spliting. This leads to a strong oscillating behavior in the event 
rates and also in the ratio, since there are slightly shifted zeroes of the oscil- 
lation term sin^(L/2-\/-B^ -|- C^). Thus there are in principle many degenerate 
points in this case as can be seen from the right hand panel. Here lines of 
constant event rates for the — transition (solid line) and for the — Vt 
transition (dashed line) in the sin^ 2^23 - e^ur plane are shown. There are two 
points were the dashed and solid hues cross. These points have exactly the 
same physical observables and are therefore not distinguishable in an exper- 
iment which uses only the total event rates. However the point in the upper 
left corner can be excluded by using the information of atmospheric neutrinos 
that sin^ 2^23 > 0.8 and e^r < 0.02 [20]. There are also many more possible 
solutions for e^,- values larger than 0.1. In order to improve the knowledge on 
e^r it is therefore necessary to include atmospheric data. 

We now come to our final results. In figures 3 and 4 we present our calculated 
NUFACT sensitivities to non-standard neutrino interactions shown as black 
solid hues for three different muon energies 20GeV, 50GeV, 100 GeV and 
150 GeV. The baseline has been chosen as 732 km. The dashed lines show the 
bounds which would be obtained by neglecting the tau mass threshold [29] 
in the cross-section. It is clearly visible that especially for low energies this 
is not a good approximation, since for example at 20 GeV one looses about 
80 % of the events. For comparison we also indicate with the grey shaded area 
the region presently excluded by the latest atmospheric data. These bounds 
are taken from [20]. The parameters were fixed as in Fig. 1: sin^ 2^23 = 0.97, 
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Figure 2. In the right hand panel hnes of constant event rates of the V/j, v^- 
transition (solid line) and of the z/^ Ur transition (dashed lines) are shown in 
the sin^ 2^23 - e^tr plane. The grey shaded band indicate the Gaussian So" errors on 
these numbers. The left hand panel shows a section across the right hand figure 
at sin^ 2^23 = 0.9. The baseline for both plots is 732 km and the muon energy is 
50GeV. All other parameters are kept fixed to: sin^ 2^13 = 0.02, sin^ 2^12 = 0.78, 
5cp = 0, Amii = 3.1 • 10-3 eV^ Amii = 3.3 • 10"^ eV^ and e|,^ = 0. 

sin2 2^i3 = 0.02, sin^ 2^12 = 0.78, 8cp = 0, Amj^ = SA-lO'^eV^, and Am^^ = 
3.3-10-5 eVl One sees that the limits on the FC parameter e^r can be improved 
by approximately two orders of magnitude by a high energy neutrino factory. 
For the low energy option (leftmost panel) the improvement in the sensitivity 
at a neutrino factory is very small compared to the present atmospheric bound. 
Note that the sensitivity on e'^^ attainable at a long baseline neutrino factory 
is worse than the present bounds by atmospheric data. 

Note also that the above bounds do not require tau charge identification. This 
is possible because the Ue-t^r transition is suppressed by sin^ 26*13 (restricted to 
be smaller than 0.1 by the Chooz experiment [36]) and ordinary matter effects 
do not come into play at the distance of 732 km considered here. In fact we 
have explicitly checked that our results are unchanged if the signs of the NSI 
parameters get reversed (in all possible combinations) with respect to what 
we have assumed. This approximation might break down at baselines longer 
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E„ = 20 GeV 




Figure 3. Here the sensitivity limits in the e^i-- e^^ plane of a neutrino factory 
compared to atmospheric neutrino data are shown for an energy resolution of 50%. 
For comparison we also indicate with the grey shaded area the region presently 
excluded by the latest atmospheric data, taken from [20]. All bounds are at 99% 
CL. The dashed line is obtained by neglecting the tau mass and is only shown for 
comparison. The black lines are calculated with the correct cross-section. All other 
parameters arc kept fixed to: values suitable to account for the present neutrino 
anomalies. Details in text. 
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Figure 4. Same as in Fig. 3 but for higher NUFACT energies. 



than 1000 km or so. 
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5 Discussion and Conclusions 

We have considered the potential of a long baseline neutrino factory in probing 
non-standard neutrino-matter interactions. We have found that the sensitiv- 
ity to flavour-changing NSI can be substantially improved with respect to 
present atmospheric neutrino data, especially at energies higher than 50 GeV 
or so, where the effect of the tau mass is small. For example, a 100 GeV NU- 
FACT can probe FC neutrino interactions at the level of few \e\ < few x 10~^ 
at 99 % C.L. The analysis we have presented requires no tau charge iden- 
tification and is based only on total event numbers, with a modest energy 
resolution at the 50 % level. In order to be useful for more refined studies a 
good detector energy resolution is required: for a 50% energy resolution the 
results are basically the same as those obtained when considering only total 
rates. It is doubtful whether a better resolution can be achieved in practice 
for the channel considered here because of hadronic tau decays. Finally note 
that the quality of the atmospheric data plays a crucial role in setting this 
limit by removing unwanted degeneracies in predicted event numbers. In con- 
trast the sensitivity on e'^^ attainable at a long baseline neutrino factory is 
worse than the present bounds by atmospheric data. The role of a NUFACT 
in probing interactions is also complementary to efforts to probe for similar 
flavour- changing effects in the charged lepton sector and has the advantage of 
being totally model independent. 
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